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Abstract 



We study the neural field equations introduced by Chossat and Faugeras in [ ] to model the 
representation and the processing of image edges and textures in the hypercolumns of the cortical area 
VI. The key entity, the structure tensor, intrinsically lives in a non-Euclidean, in effect hyperbolic, 
space. Its spatio-temporal behaviour is governed by nonlinear integro-differential equations defined 
on the Poincare disc model of the two-dimensional hyperbolic space. Using methods from the theory 
of functional analysis we show the existence and uniqueness of a solution of these equations. In 
the case of stationary, i.e. time independent, solutions we perform a stability analysis which yields 
important results on their behavior. We also present an original study, based on non-Euclidean, 
hyperbolic, analysis, of a spatially localised bump solution in a limiting case. We illustrate our 
theoretical results with numerical simulations. 
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O 1 Introduction 

Chossat and Faugeras in [ ] have introduced a new and elegant approach to model the processing of 
image edges and textures in the hypercolumns of area VI that is based on a nonlinear representation 
of the image first order derivatives called the structure tensor. They assumed that this structure tensor 
was represented by neuronal populations in the hypercolumns of VI that can be described by equation 
similar to those proposed by Wilson and Cowan [26]. 

Our investigations are motivated by the work of Bressloff, Cowan, Golubitsky, Thomas and Wiener [8, 9 
on the spontaneous occurence of hallucinatory patterns under the influence of psychotropic drugs and 
the further studies of Bressloff and Cowan [7, (i, ■">]. We hardly think that the natural spatial extension 
of our model would lead to an exciting anlysis of hyperbolic hallucinatory patterns. But, this requires 
first to better understand the a-spatial model and this is the subject of this present work. The a-spatial 
model can also be linked to the work by Ben-Yishai [3] and Hansel, Sompolinsky [ ] on the ring model 
of orientation. 

The aim of this paper is to present a general and rigorous mathematical framework for the modeling of 
neuronal populations in one hypercolumn of VI by the structure tensor which is based on miscellaneous 
tools of functional analysis. We illustrate our results with numerical experiments. In section 2 we briefly 
introduce the equations, in section 3 we analyse the problem of the existence and uniqueness of their 
solutions. In section 4 we deal with stationary solutions. In section 5, we present an analysis of what 
we called a hyperbolic radially symmetric stationary-pulse in a limiting case. In the penultimate, we 
present some numerical simulations of the solutions. We conclude in section 7. 



2 The model 

We recall that the structure tensor is a way of representing the edges and textures of a 2D image T(x, y) 
[4, 21]. Moreover, a structure tensor can be seen as a 2 x 2 symmetric positive matrix. 
We assume that a hypercolumn of VI can represent the structure tensor in the receptive field of its 
neurons as the average membrane potential values of some of its membrane pouplations. Let T be a 
structure tensor. The average potential V(T, t) of the column has its time evolution that is governed 
by the following neural mass equation adapted from [ ] where we allow the connectivity function W 
to depend upon the time variable t and we integrate over the set of 2 x 2 symmetric definite-positive 
matrices: 

d t V(T,t) = -aV(T,t) + J sm{2) W(T,T',t)S(V(T',t))dr + I(T,t) Vt > 

V(T,0) = V Q (T) { ' 

The nonlinearity S is a sigmoidal function which may be expressed as: 

S(x) l 



1 + e-^ x 



where fj, describes the stiffness of the sigmoid. / is an external input. 

The set SPD(2) is the set of 2 x 2 symmetric positive-definite matrices. It can be seen as a foliated 

manifold by way of the set of special symmetric positive definite matrices SSPD(2) = SPD(2) n SL(2, R). 

Indeed, we have: SPD(2) h ™ SSPD(2) x R+. Furthermore, SSPD(2) i = ? " D, where D is the Poincare 

Disk, see e.g. [11]. As a consequence we use the following foliation of SPD(2): SPD(2) = D x R+, 
which allows us to write for all T G SPD(2), T = (z, A) with (z, A) € D x M+. T, z and A are related by 
the relation det(T) = A 2 and the fact that z is the representation in D of T € SSPD(2) with T = AT. 
It is well-known [20] that D (and hence SSPD(2)) is a two-dimensional Riemannian space of constant 
sectional curvature equal to -1 for the distance noted d.2 defined by 

\z - z'\ 

d,2(z,z') = arctanh- — -. 

1 1 — zz I 

It follows, e.g. [24, 11], that SDP(2) is a three-dimensional Riemannian space of constant sectional 
curvature equal to -1 for the distance noted d defined by 



d (T,T) - ^2(logA-logA') 2 + rfi(z,z') 



As shown in proposition (A. 0.1) of appendix A it is possible to express the volume element dT in 
(zi, 02) A) coordinates with z = z\ + iz^: 



«■-»*%<£$? 



We rewrite (1) in (z, A) coordinates: 



d t V(z, A, t) = -aV(z, A,t) + 8V2 / W(z, A, z', A', t)S(V(z', A', t))— ) 2 + I(z, A, t) 

Jo Jd a U — \ z I ) 



l ' +ca r dA' fly' rly' 

We get rid of the constant 8V2 by redefining W as 8\/2W. 

d t V(z, A, t) = -aV(z, A, t) + J Q + °° J B W(z, A, z', A', t)S(V(z>, A',t))^- (1 ^ffi 2)2 + Hz, A, t) 

Vt>0 
V(z,A,0) = V (z,A) 

(2) 
The aim of the following sections is to establish that (2) is well-defined and to give necessary and sufficient 
conditions on the different parameters in order to prove some results on the existence and uniqueness of 
a solution of (2). 



3 The existence and uniqueness of a solution 

The aim of this section is to give theoretical and general results of existence and uniqueness of a solution 
of (1). In the first subsection 3.1 we study the simpler case of homogeneous solutions of (1), i.e. of 
solutions that are constant with respect to the tensor variable T. We then we study in 3.2 the useful 
case of the semi-homogeneous solutions of (1), i.e. of solutions that depend on the tensor variable but 
only through its z coordinate in D, and we end up in 3.3 with the general case. 

3.1 Homogeneous solutions 

A homogeneous solution to (1) is a solution V that does not depend upon the tensor variable T for a 
given homogenous input I(t) and a constant initial condition Vq. In (z, A) coordinates, a homogeneous 
solution of (2) is defined by: 

V(t) = -aV(t) + W(z, A, t)S(V(t)) + I(t) 

where: 

^• A ' t) = r7„ H ' (2 ' ay ' a ' ,) f(i f q#? (3) 

Hence necessary conditions for the existence of a homogeneous solution are that: 

• the double integral (3) is convergent, 

• W(z,A,t) does not depend upon the variable (z,A). In that case, we note W(t) instead of 
W(z,A,t). 

In the special case where W(z,A,z',A',t) is a function of only the distance do between (z, A) and 
(*',A'): 

W(z, A, z', A', t) = w (j2(logA-logA')*+<q(z, z '), t\ 

the second condition is satisfied. We postpone the verification of this fact until the following section. To 
summarize, the homogeneous solutions satisfy the differential equation: 

V(t) = -aV{t)+W(t)S(V(t))+I(t) t>0 ,, 

V(0) = V [ ' 

3.1.1 A first existence and uniqueness result 

Equation (2) defines a Cauchy problem and we have the following theorem. 

Theorem 3.1.1. If the external current 1 (t) and the connectivity function W(t) are continuous on some 
closed interval J containing 0, then for all Vq in R, there exists a unique solution of (4) defined on a 
subinterval Jo of J containing such that V(0) = Vq. 

Proof. It is a direct application of Cauchy' theorem on differential equations. We consider the mapping 
f : J xR^R defined by: 

f(t, x) = -ax + W(t)S(x) + I(t) 

It is clear that / is continuous from J x R to I. We have for all x,y GR and t E J: 

\f(t,x) - f(t,y)\ < a\x -y\ + \W(t)\S'Jx - y\ 

where S^ = SUp xeR \S'(x)\. 

Since, W is continuous on the compact interval J, it is bounded there by C > and: 



\f(t,x)-f(t,y)\<(a + CS' m )\x-y\ 



We can extend this result to the whole time real line if / and W are continuous on 



□ 



Proposition 3.1.1. If I and W are continuous on M., then for all Vq in M., there exists a unique solution 
of (4) defined on M. such that V(0) — Vq. 

Proof. We have already shown the following inequality: 

\f(t,x)-f(t,y)\<a\x-y\ + \W(t)\S'Jx-y\ 

Then / is locally Lipschitz with respect to its second argument. Let V be a maximal solution on Jo and 
we denote by /3 the upper bound of Jo- We suppose that /3 < +00. Then we have for all t > 0: 

V(t) = e- at V + f e- a(t - u) W(u)S(V(u))du+ f e- a(t - u) I{u)du 
Jo Jo 

^\V(t)\<\V \+S m e au \W(u)\du+ e au \I(u)\du Wg[0,/3] 
Jo Jo 

where S m = sup.,, GR |5(a;)|. 

But theorem B.0.2 ensures that it is impossible, then /3 = +00. The same proof with the lower bound of 

Jo gives the conclusion. □ 

3.1.2 Simplification of (3) in a special case 

Invariance In the previous section, we have stated that in the special case where W was a function of 
the distance between two points in D x Rf, then W(z,A,t) did not depend upon the variables (z, A). 
We now prove this assumption. 

Lemma 3.1.1. When W is only a function of do(T,T'), then W does not depend upon the variable T. 

Proof. We work in (z, A) coordinates and we begin by rewritting the double integral (3) for all (z, A) G 



x 



W(z,A,t)=f J W (y/2(log A -log A')* + 4(z,z'),t\ 



The change of variable A' — > t yields 



A 

W(z,A,t) = J JjvU2{\ogA>y + dl{z,z>),t\ 



dA' dz[dz' 2 

^(l-lz'l 2 ) 2 

dA' dz[dz' 2 
1^ {l-\z'\ 2 ) 2 



And it establishes that W does not depend upon the variable A. To finish the proof, we show that the 
following integral does not depend upon the variable z G D: 

^)=/ D «<M^)>(I^ <5> 

where / is a real-valued function such that S(z) is well defined. 

We express z in horocyclic coordinates: z — n s a r .O (see appendix D) and (5) becomes: 

E(z) = / / f(d2(n s a r .O,n s ia r '.0))e~~ 2r ds'dr' 



f{d<2,{n s - s ia r .O,a r i.O))e 2r ds'dr' 
With the change of variable s — s' = —xe 2r , this becomes: 

3{z)= [ [ f(d 2 (n_ xe 2 r a r .O,a rl .O))e- 2{r '- r) dxdr' 

JrJm 

The relation n_ xe 2 r a r .O = a r n~ x .O (proved e.g. in [ ]) yields: 

HO) = / J f{d 2 {a r n^ x .O,a r ,.0))e- 2 ^- r ^'dxdr' 



f{d 2 (0,n x a r ,- r .O))e- 2< - r '- r) dxdr' 
f(d 2 (0,n x a y .O))e~ 2v dxdy 



= / f(d 2 (0,z'))dm(z') 

Jb 

lz' 2 



dz' dz' 
with z' = z[ + iz' 2 and dm(z') = : — r^r^, which shows that S(z) does not depend upon the variable 



z. 



as announced. □ 



Mexican hat connectivity In this paragraph, we push further the computation of W in the special 
case where W does not depend upon the time variable t and takes the special form suggested by Amari in 
[ ], commonly referred to as the "Mexican hat" connectivity. It features center excitation and surround 
inhibition which is an effective model for a mixed population of interacting inhibitory and excitatory 
neurons with typical cortical connections. It is also only a function of e?o(T, T'). 
In detail, we have: 



W(z, A, z'A') = w I y'2(log A - log A') 2 + <%{z, z') 

where: 

x 2 x 2 

i \ * ^a 2 A ~~a 2 

w(x) = e °i e °2 

v/2^2 y^TTof 

with < (Ti < cr 2 and < A < 1. 

In this case we can obtain a very simple closed-form formula for W as shown in the following lemma. 

Lemma 3.1.2. When W is a mexican hat function of d 2 (T 1 T') and independent oft, then: 

3 

W=7 Y ( CT i e2CT W(^ / 2~ (7 i) - Aa 2 e 2 ^erf(y2a 2 ^ (6) 

where erf is the error function defined as: 

erf(x) = —= \ e~ u du 
V 71- Jo 

Proof. The proof is given in appendix C. □ 

3.2 Semi-homogeneous solutions 

A semi-homogeneous solution of (2) is defined as a solution which does not depend upon the variable A. 
In other words, the populations of neurons are not sensitive to the determinant of the structure tensor. 
The neural mass equation is then equivalent to the neural mass equation for tensors of unit determinant. 
We point out that semi-homogeneous solutions were previously introduced by Chossat and Faugeras in 
[11]. They also performed a bifurcation analysis of what they called H-planforms. In this section, we 
define the framework in which their equations make sense. 

d t V(z,t) = -aV(z,t) + j W sh (z,z',t)S(V(z',t))dm(z')+I{z,t) t>0 

V(z,0) = V (z) [ '' 

where 

r*+oo 



W sll {z,z',t) = 



r + °° dA' 

J W{z,A,z\A\t) a ~ r 



We have implicitly made the assumption, that W sh does not depend on the coordinate A. Some condi- 
tions under which this assumption is satisfied are described below. 

We now deal with the problem of the existence and uniqueness of a solution to (7) for a given initial 
condition. We first introduce the framework in which this equation makes sense. 



3.2.1 The well-posedness of equation (7) 

Let J be an open interval of R. We assume that: 

• (CI): V(z, z', t) G D 2 x J, W sh (z, z', t) = w sh {d 2 (z, z'),t), 

• (C2): W sh G C (J,£ X (B)) where W sh is defined as W sh (z,t) = w sh (d 2 (z,0),t) for all 

(z,t) gDx J, 

. (C3): sup feJ ||W s/l (t)|| L1 < +cx> where ||W s/l (t)|| L1 = f J B \W sh (d 2 (z,0), i)|dm(z). 

Note that conditions (C1)-(C2) and lemma 3.1.1 imply that for all z G D, / D |W s,l (z, z',i)| dm(z') = 
l|W s/l (t)|| L i. And then, for all z G O, the mapping z' — ► W sh (z, z',t) is integrable on D. We deduce 
that for all (z,t) eDx J: 

< f \W sh (z ) z',t)\\S(V(z',t))\dM^) 



W sh (z,z' ',t)S(V(z' ,i))dm(z') 

<S m [ \W sh (z,z',t)\dm(z') <5 m sup||W sft (t)|| i i <+oo 
jd te.j 

The last inequality is a consequence of lemma 3.1.1 which shows that J D |W s,I (z, z', t)\ dm(z') is not a 
function of z. 

Finally for all z £ D and t £ J, the righthand side of equation (7) is well-defined. 
We introduce the following mapping, defined on J x J 7 , where J 7 is a yet to be defined functional space: 

/ : (t, (/>) -> /(t, 0) such that /(*, 0)(z) = / M^ s,l (z, z', t)S'(0(z'))dm(z') (8) 



Our aim is to find a functional space J- where (8) is well-defined and the function / maps J- to J- for all 

is. 

A natural choice would be to choose as a L p (B,<im)-integrable function of the space variable with 

1 < p < +oo. Unfortunately, the homogeneous solutions (constant with respect to z) do not belong 

to that space. Moreover, a valid model of neural networks should only produce bounded membrane 

potentials. That is why we focus our choice on functional spaces of bounded, or essentially bounded, 

functions. 

The well-posedness of equation (7) in the case T = L°°(D) Our first choice is T = L°°(D). 

The Fischer-Riesz's theorem ensures that L°°(D) is a Banach space for the norm: || VII oo = inf{C > 
0, |V , (^)| < C for almost every z G D}. We have the following proposition. 

Proposition 3.2.1. IJW sh satisfies conditions (C1)-(C3) then f is well defined and is from JxL°°(B) 
to L°°(0). 

Proof. For all <f> G T and t G M, fit, <fi) is obviously measurable and we have from (8): 

||/(* > 0)||oo<S m sup||W h (t)|Ui<+oo 

and hence f(t, <p) G L°°(D). □ 

The well-posedness of equation (7) in the case T = C&(B) We now choose J- — Cb(D) the space 
of the bounded continuous functions on D. As shown in, e.g. [18], this is a Banach space with respect 
to the uniform norm: HV'llcb(B) — su Pzgb IV'C 2 )!- The previous proposition 3.3.1 still holds: 

Proposition 3.2.2. IfW sh satisfies conditions (C1)-(C3) then f is well defined and is from J x C&(B) 
toC b (p). 



Proof. For all 4> £ C b (D) and t £ R: 

ll/(t,0)lk(D) < ^supHW^COIUi < +^> 
teJ 

so /(£, 0) is bounded. It remains to show that for all (t,z,4>) €E J x D x C&(0) the mapping z -4 /(£, </>)(z) 
is continuous on D. We fix t € J and write z in horocyclic coordinates z — n s a r .O with (s,r) £ R 2 and 
define / on R 2 as: 

/: (s,r)^ f [ W sh (d 2 {n s a r .O,n s ,a r , .0),t)S((t>(n s ,a r , .0))e- 2r ' ds'dr' 

Jr Jm 

A method similar to the one used in section 3.1.2 leads to: 

/(s,r)= / [ W sh (d 2 (O,n x a y .O),t)S(0(n s+xe 2ra r+y .O))e~ 2y dxdy 

JK.JM 

We now use a theorem on integrals depending on a parameter. It is easy to verify that 

1. for all (s,r) £ R 2 the function (x,y) — > W sh (d 2 (0,n x a y .O),t)S(4>(n s+xe 2ra r+y .O))e~ 2y is measur- 
able on R 2 , 

2. for almost every (x,y) £ R 2 the function (s,r) — > W sh (d 2 (0,n x a y .O),t)S(4>(n s+xe 2 r a r+y .O))e~ 2v 
is continuous on R 2 , 

3. for all (s,r) £R 2 , 

\W sh (d 2 (0,n x a y .O),t)S(cj ) (n s+xe 2 r a r+y .O))e- 2 y\<\W sh (d 2 (0,n x a y .O),t)\S m 

and (x, y) — > \W sh [d 2 (0,n x a y .O),t)\ is integrable on R 2 . 

It follows that the function / is continuous on R 2 and / is continuous on D. 

Finally /(£, <j>) belongs to C b (D). D 

3.2.2 The existence and uniqueness of a solution of (7) 

From now on, J 7 is a functional Banach space for the norm || • \\jr. We suppose that all the hypotheses are 
verified so that / is well-defined from J x T to J- with J an open interval containing 0. In the previous 
section we have already presented two different examples for T: L°°(D) and C(,(D). 
We rewrite (7) as a Cauchy problem defined on T: 

V'(t) = -aV(t) + f(t,V(t)) + I(t) t£j 

no) = v w 

Theorem 3.2.1. If the external current I belongs to 0^,7^) with J an open interval containing and 
W sh satisfies conditions (C1)-(C3), then for all Vq £ J- , there exists a unique solution of (9) defined on 
a subinterval Jq of J containing 0. 

Proof. We prove that / is continuous on J x T . We have 

f{t, x) - f(s, y) = f W sh {-,z', t) (S(x) - S(y))dm(z') 

Jo 

+ J (w sh (-,z', t) - W s % z', s^j S(y)dm(z'), 

\\f(t,x) - f(s,y)\\r < S' m W \\x -y\\ F + S m \\W s \t) - W s/l ( s )|| L1 

Because of condition (C2) we can choose \t — s\ small enough so that 

II ~W sh (t) — ~W sh (s) || . 1 is arbitrarily small. This proves the continuity of /. Moreover it follows from the 

previous inequality that: 

\\f(t,x)-f(t,y)\\r<S' m W \\x-y\\r 

with Wo = sup te j ||W s,l (t)|| i i. This ensures the Lipschitz continuity of / with respect to its second 
argument, uniformly with respect to the first. The Cauchy-Lipschitz theorem on a Banach space yields 
the conclusion. □ 



This solution, denned on the subinterval J of R can in fact be extended to the whole real line, and 
we have the following proposition. 

Proposition 3.2.3. If the external current I belongs toC(M.,J 7 ) andW sh satisfies conditions (C1)-(C3) 
with J = R, then for all Vq G J- ' , there exists a unique solution of (9) defined on R. 

Proof. The proof is a direct application of theorem B.0.3. If / is Lipschitz continuous with respect to its 
second argument, the righthand side of (9) is also Lipschitz continuous. Now / is Lipschitz continuous 
because 

\\f(t,x) - f(t,y)\\? < S' m W \\x -y\\r Vi G R 

D 

3.2.3 The boundedness of a solution of (7) 

We assume that J 7 is a Banach space chosen so that the mapping / is well defined from K x J to J. 
Then the following propostion holds. 

Proposition 3.2.4. If the external current I belongs to C(R + , J-) and is bounded in time, i.e. sup tgR + ||/(£)||jf < 
+oo, and W sh satisfies conditions (C1)-(C3) with J = R + * , then the solution of (9) is bounded for each 
initial condition Vq G T . 

Proof. For all t G R + * we integrate (7) over [0,t]: 

V(z,t) = e~ at V (z)+ [ e - a{t - u) [ W sh (z,z' ,u)S(V(z' ,u))dm(z')du+ [ e- a{t - u) I(z,u)du 
Jo Jo Ja 

The following upperbound holds 

\\V(t)\\r < ^"l^oll^ + - (s m W + sup \\I(t)\\A (1 - e- at ) , (10) 

and hence 

\\V(t)\\r < \\V \\jr + - (s m W + sup ||7(t)||^ 

a \ teM+ 

which shows that the solution is bounded for each initial condition Vb G T . 

U 

The upperbound (10) yields a simple attracting set for the dynamics of (7) as shown in the following 
proposition. 

Proposition 3.2.5. Let p = £ (S ,m Wo + sup teR + ||I(i)||^r). The open ball B p of T of center and 
radius p is stable under the dynamics of equation (7). Moreover it is an attracting set for this dynamics 
and ifVotf: B p and T = inf {t > such that V(t) G B p } then: 

T<ilog^"- 



a \ p 

Proof. We can rewrite (10) as: 

\\V(t)y < e- at (\\V y - - (s m W + sup ||/(t)||^) + - (s m W + sup ||I(t)||^ 



'•{m\r-l) + l (ID 



If V Q G B p this implies ||V(t)||jr < f (1 + er at ) for alii > and hence ||F(i)||^ < p for all i > 0, proving 
that B p is stable. Now assume that ||V(i)||^- > p for all i > 0. The inequality (11) shows that for i large 
enough this yields a contradiction. Therefore there exists io > such that ||^(io)||jr = P- At this time 
instant we have 

and hence 

i0 <I log W,||.-p 



a \ p 

n 



3.3 General solution 

We now deal with the general solutions of equation (1). We first give some hypotheses that the con- 
nectivity function W must satisfy. We present them in two ways, first on the set of structure tensors 
considered as the set SPD(2) and, second on the set of tensors seen as D x K.+ . Let J be a subinterval 
of E. We assume that: 

• (HI): V(T,T',t) € SPD(2) x SPD(2) x J, W(T,T',t) = W(d (T,V),t), 

• (H2): W e C (J,L 1 (SPD(2))) where W is defined as W(T,t) = W (do(T,Id 2 ),t) for all (T,t) G 
SPD(2) x J where Id 2 is the identity matrix of M 2 (M), 

. (H3): Vi E J, sup teJ ||W(t)||ii < +00 where ||W(i)|| L1 */ / gpD(2) |W(d (T,Id 2 ),i)|dT. 
We now express these hypotheses for the representation in (z, A) € D x E+ of structure tensors: 

• (Hlbis): V(z, z', A, A', i) e D 2 x (E+) 2 x E, W(z, A, z', A', i) = W(da(z, «0> I log(A) - log(A')|, t), 

. (H2bis): W € C (J,^ 1 (D x E+)) where W is defined as W(z, A,i) = W (d 2 (z,0),\log(A)\,t) for 

allV(z,A,i) € DxE+ x J, 

• (H3bis): Vi € J, sup teJ ||W(i)|| L i < +00 where 

\\W(t)\\»= f I \W(d 2 (z,0),\\og(A)\,t)\^dm(z) 



3.3.1 Functional space setting 

We need to settle on the choice of a Banach functional space JF for the membrane potential as in section 
3.2. Our study of the semi-homogeneous case suggests the following choice: T = L°°(D x E+). As 
D x E+ is an open set of E 3 , J 7 is a Banach space for the norm: ||0||jf- = sup zSD sup AeM + \4>{z. A)|. 
We introduce the following mapping f 9 : (t, <fi) —> f 9 (t, <j>) such that: 



/»(*, 4>){z, A) = -a<t>(z, A)+ / W( d 2 (z, z') 



!°g(^) 



, t\ S (<j>(z', A')) ^dm(z')+/( Z , A, t) (12) 



Proposition 3.3.1. If I E C(J,J-^J with sup teJ ||J(t)||jr < +00 and W satisfies hypotheses (Hlbis)- 
(H3bis) then f 9 is well-defined and is from J x T to J- '. 



Proof. V(z, A, t) £ D x E+ x E, we have: 



W[d 2 {z,z'), 



tog(^) 



dA' 



t 5(^',A'))-dm( Z ') 



A' 



<S m sup||W(t)|| i i <+oo 
te.J 



D 



3.3.2 The existence and uniqueness of a solution of (2) 

We rewrite (2) as a Cauchy problem: 



f d t V(z 7 A,t) = -aV(z,A,t) + J BxRt W(d 2 (z,z'),\\og(^)\,t)S(V(z\A',t)) d ^dm(z')+I(z,A,t) 

I V(z,A,0) = V (z ) A) 

(13) 

Theorem 3.3.1. If the external current I belongs to C(J,J-) with J an open interval containing and 
W satisfies hypotheses (Hlbis)- (H3bis), then fo all Vq G J- ', there exists a unique solution of (13) defined 
on a subinterval Jq of J containing such that V(z, A, 0) = Vq(z, A) for all (z, A) € D x E+. 



Proof. The proof is a direct adaptation of the proof of theorem 3.2.1. □ 

Proposition 3.3.2. If the external current I belongs to C(K, J-) and W satisfies hypotheses (Hlbis)- 
(HSbis) with J = M., then for all Vq £ J- , there exists a unique solution of (13) defined on M. such that 
V(z, A, 0) = V (z, A) for all (z, A) <= D x R+ . 

Proof. The proof is readily adapted from that of proposition 3.2.3. □ 

3.3.3 The intrinsic boundedness of a solution of (2) 

In the same way as in the homogeneous case, we show a result on the boundedness of a solution of (2). 
The proofs of the following properties are exactely the same as those in section 3.2.3. 

Proposition 3.3.3. If the external current I belongs to C(R + , J-) and is bounded in time sup tGR + ||/(£)||jf < 
+oo and W satisfies hypotheses (Hlbis)-(HSbis) with J = R + , then the solution of (13) is bounded for 
each initial condition Vq G J-. 

If we set: 

p9 d g l { S m W$ + sup \\I(t)\\r) 

a teR+ 

where Wg = sup teR + ||W(t)||^i. The following corollary is a consequence of the previous proposition. 
Corollary 3.3.1. IfV £ B p9 and T CJ = inf{£ > such that V{t) e B pg } then: 

T 9 <llogf 2 "^-^ 
a \ p 9 

4 Stationary solutions 

We look at the equilibrium states, noted V® of (2), when the external input I and the connectivity W 
do not depend upon the time. We assume that W satisfies hypotheses (Hlbis)-(H2bis). We redefine for 
convenience the sigmoidal function to be: 

S(x) ' 



1 + e~ x ' 
so that a stationary solution (independent of time) satisfies: 



= -aV°(z, A) + f W (d 2 (z, z'), 

Jo V 



ks(£) 



rIA' 

S(nV°{z', A'))—dm{z') + I(z, A) (14) 



'A 
We define the nonlinear operator from JF to J 7 , noted Q^, by: 

G,(V)(z, A) = / W (d 2 (z, z'), log(^) ) SQiVtf, A'))^-dm(z') (15) 

Finally, (14) is equivalent to: 

aV°(z,A) = g^V)(z,A) + I(z,A) 

4.1 Study of the nonlinear operator Q^ 

We recall that we have set for the Banach space T = L°°(D x R+) and proposition 3.3.1 shows that 
Q^ . : T — > J ' . We have the further properties: 

Proposition 4.1.1. Q^ satisfies the following properties: 

• WGpiVi) - Gv(V 2 )\\r < nWgSlnWVi - VaWr for all /i > 0, 

• jJ> —> Qjj, is continuous on R + , 
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Proof. The first property was shown to be true in the proof of theorem 3.2.1. The second property 
follows from the following inequality: 



IISmi (Y) - 9^{V)\\ r < |/ii - ^\W°S' m \\V\ 



T 



D 

We denote by Gi and Goo the two operators from J 7 to J 7 defined as follows for all V £ T and all 

(z,A) eBxR+: 



Gi(V)(z,A)= [ w(d 2 (z,z'), 

JDxRt \ 



!og(^) 



dA' 

V{z , ,A')^dm{z') ) 



(16) 



and 



0oo(V)(«,A) 



W(d 2 {z,z'), 



log(^) 



)if(F(^,A'))^dm(^) 



where H is the Heaviside function. 

It is straightforward to show that both operators arc well-defined on T and map T to T ' . Moreover 
the following proposition holds. 

Proposition 4.1.2. We have 

Ga > Goo 

Proof. It is a direct application of the dominated convergence theorem using the fact that: 

S(fiy) — ► H(y) p.p y£R 

□ 



4.2 The convolution form of the operator Q^ in the semi-homogeneous case 

It is convenient to consider the functional space J- sh = L°°(D) to discuss semi-homogeneous solutions. 
A semi-homogeneous persistent state of (2) is deduced from (14) and satisfies: 



where the nonlinear operator Gt h from F sh to F sh is defined for all V £ J- sh and z £ D by: 

Gf(V)(z) = f W sh (d 2 (z,z'))S(nV(z'))dm(z') 



(17) 



We define the associated operators, Gf h ,Gf h ' 



sh (~ir\f \ / jxrsh 



0p(V)(*) = / W sh (d 2 (z,z'))V(z')dm(z') 



G^(V)(z)= / W sh (d 2 (z,z'))H(V(z / ))dm(z / ) 



We rewrite the operator Gt in a convenient form by using the convolution in the hyperbolic disk. First, 
we define the convolution in a such space. Let dg denote the Haar measure on the group G = SU(1, 1) 
(see [17]), normalized by: 

def 



f(g-0)dg^ / f(z)dm(z), 
for all functions of L : (D). Given two functions f%, f 2 in L 1 (D) we define the convolution * by: 

(A*/ 2 )(z)= I f 1 (g-0)f 2 (g- 1 -z)dg 
Jg 

We recall the notation W sh (z) d = W sh (d 2 (z, O)). 
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Proposition 4.2.1. For all /z > and V 6 J rs ' 1 we have: 

gf(V) = W h * S{fiV), Gf h (V) =W sh *V and G 8 £{V) = r' 1 * H(V) (18) 

Proof. We only prove the result for Q^. Let zfD, then: 



0* h (V)(*) = /_W sh (d 2 (^^)),S(MF(zO)dm(/) = \w sh {d 2 (z,g-0))S^V{g-0))dg 

W sh (d 2 (gg' 1 -z,g- 0))S{^V{g ■ 0))dg 



IG 

and for all g G SU(1,1), d-i{z,z') — d 2 (g ■ z,g ■ z') so that: 

G?{V){z) = f W sh (d 2 (g- 1 ■ z, 0))S(fiV(g ■ 0))dg = W sh * S(jjV){z) 
Jg 



D 



Let b be a point on the circle SB. For z £ D, we define the "inner product" < z, b > to be the 
algebraic distance to the origin of the (unique) horocycle based at b through z (sec [ ]). Note that 
< z,b > does not depend on the position of z on the horocycle. The Fourier transform in D is defined 
as (see [17]): 

h(X, b) = / h{z)e { ~ lX+1)<z ' b> dm(z) V(A, b) G R x <9D 



for a function ft, : D — ► C such that this integral is well-defined. 

Lemma 4.2.1. T/ie Fourier transform in D, VF (A, 6) o/ W 8 does not depend upon the variable b G <9B. 
Proof. For all A € R and 6 = e i0 G <9B, 

W S ''(A,6) = / W sh (z)e { - lX+1}<z ' b> dm{z). 

Jb 

We recall that for all £ R r^ is the rotation of angle <\> and we have W s (rj, ■ z) = W s (z), dm(z) — 
dm^cf, • z) and < z, b >=< r^ ■ z, r^ ■ b >, then: 

W h (X,b) = / W s/l (r_ ■ z)eS- lX+ ^ <r - e - z ' 1> dm(z) = / W s ' i (z)e^ a+1 ) <2 ' 1> dm(z) = J W'V) 

D 

We now introduce two functions that enjoy some nice properties with respect to the Hyperbolic 
Fourier transform and are eigenfunctions of the linear operator Qf . 

Proposition 4.2.2. Let e x , b {z) = e (-<A+i)<*,&> and $ A (z) = J™ e^ x+ ^ <z ' b> db then: 



• G( h (ex, b ) = W h (\)e x , b 

• g? h (®x) = w sh (\)$ x 

Proof. We begin with 6 = 1 G 9D and use the horocyclic coordinates. We use the same changes of 
variables as in lemma 3.1.1: 

gj h (ex,i)(n 3 at.O) = f W sh (d 2 (n s a t ■ 0,n s ,a t , ■ 0))e^ iX -^' dt'ds' 

JR 2 

W sh {d 2 {n s - S ,a t ■ 0,Of • 0))e ( - lX ~ l)t ' dt' ds' 



W sh (d 2 {a t n^ x ■ O, a v ■ 0))e ( ~ lX - x)t ' + 2t dt! ' dx = / W sh (d 2 (0, n x a V -t ■ O)^- 1 ^ 1 ' +2t dt' dx 

Jr 2 
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W sh (d 2 (0,n x a y ■ 0))e(- lX ^y+V +2t dydx = e (-*A+l)<».» t .O > l> W "' , ( A ) 



By rotation, we obtain the property for all b £ 90. 

For the second property [17, Lemma 4.7] shows that: 



W sh * $ A 



(z) = f e(' iA+1 > <z ' 6> W Sh (A)d& = $ A (z)W S "(A) 

JdO 



D 



A consequence of this proposition is the following lemma. 



Lemma 4.2.2. The linear operator Qf h is not compact and for all \i > 0, the nonlinear operator Qf/ 1 is 
not compact. 

Proof. The previous proposition 4.2.2 shows that Qf has a continuous spectrum which iimplies that is 
not a compact operator. 

Let U be in F sh , for all V £ J- sh we differentiate Off 1 and compute its Frechet derivative: 

D (gf) v (V)(z) = / W° h (d 2 (z, z'))S'(U(z'))V(z')dm(z') 

Jo 

If we assume further that U does not depend upon the space variable z, U(z) = U$ we obtain: 

D(g?) Uo (v)(z) = s'(u )gf h (v)(z) 

If Q^ 1 was a compact operator then its Frechet derivative D [p^ h j would also be a compact operator, 
but it is impossible. As a consequence, Q s Jf is not a compact operator. □ 

4.3 The convolution form of the operator Q^ in the general case 

We adapt the ideas presented in the previous section in order to deal with the general case. We recall 
that if H is the group of positive real numbers with multiplication as operation, then the Haar measure 

dx 
dh is given by — . For two functions fi, f 2 in L 1 (D x R+) we define the convolution * by: 
x 

(A * f 2 )(z, A) *! [ [ fi(g-0,h- lJ/aGT 1 ■ z, h' 1 ■ A)dgdh 

We recall that we have set by definition: W(z, A) = W(d 2 (z, 0), | log(A)|). 
Proposition 4.3.1. For all fi > and V £ T we have: 

g„{V) = W*S(/iV), Gi(V) =W*V andg^V) = W*H(V) (19) 

Proof. Let (z, A) be in D x R+. We follow the same ideas as in proposition 4.2.1 and prove only the first 
result. We have 

d/\' 
g,(V)(z,A) = J 0x ^Wfa(z,z'),\] O g&)\)S(jMV(z',A'))? Er dmW 

dA' 

= Ig kt W (^GT 1 • z, O), | log (A) |) S (vV(g ■ O, A')) dg — 

= IgIh w (^Gr 1 ■ z,0),\log(h- 1 ■ A) |) S{nV{g-0,h-l))dgdh 
= W*S{nV)(z,A) 



D 



We next assume further that the function W is separable in z, A and more precisely that W(z, A) 
Wi(«)W 2 (bg(A)) whereWjfz) = W x {d 2 {z,Q)) and W 2 (log(A)) = W 2 (|log(A)|) for all (2, A) £ DxK; 
The following proposition is an echo to proposition 4.2.2. 
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Proposition 4.3.2. Let e x , b {z) = e (- iA + 1 )< z - b >, $ A ( Z ) = i / SD e( lA+1 ) <2 ' 6> d& and h^(A) = e ^ log ( A ) 
i/ien; 

• SKeA,b^)=Wi(A)W 2 (e)e A ,b^ 

• ft(*A^)=Wi(A)W 2 (0*A^ 

where W2 is i/ie usual Fourier transform of W2 . 

Proof. The proof of this proposition is exactly the same as for proposition 4.2.2. Indeed: 

, rlA' 
a,(e A , 6 fc*)(*,A) = W 1 * eA , 6 (z)/ R +W 2 (log(A))e^^(A)_ 

= Wi * e x , b (z) (J M W 2 {y)e-^ydy) e***^ 

D 

A straightforward consequence of this proposition is an extension of lemma 4.2.2 to the general case: 

Lemma 4.3.1. The linear operator Qf h is not compact and for all /j, > 0, the nonlinear operator Q* h is 
not compact. 

4.4 The set of the solutions of (14) 

Let B^ the set of the solutions of (14) for a given slope parameter fj,: 

B^ = {VeT\-aV + 0„(V) + / = 0} 
We have the following proposition. 

Proposition 4.4.1. If the input current I is equal to a constant Iq, i.e. does not depend upon the 
variables (z, A) then for all \x G M + , B^ 7^ 0. In the general case I G J- , if the condition ^,S' m W^ < a is 
satisfied, then Card(B^) = 1. 

Proof. Due to the properties of the sigmoid function, there always exists a constant solution in the case 
where I is constant. In the general case where I G T , the statement is a direct application of the Banach 
fixed point theorem, as in [13]. □ 

Remark 4.4.1. It should be clear that if the input current does not depend upon the variables (z, A) and 
if the condition (j,S' m WQ < a is satisfied, then there exists a unique stationary solution which, in effect, 
does not depend upon the variables (2, A). 

If on the other hand the input current does depend upon these variables, is invariant under the action 
of a subgroup of U (1,1), the group of the isometries o/D (see (D)), and the condition fiS' m WQ < a is 
satisfied, then the unique stationary solution will also be invariant under the action of the same subgroup. 
When the condition [xS^Wq < a is satisfied we call primary stationary solution the unique solution in 
B M . 

4.5 Stability of the primary stationary solution 

In this subsection we show that the condition [iS^Wq < a guarantees the stability of the primary 
stationary solution to (2). 

Theorem 4.5.1. We suppose that I G J- and that the condition ^,S' m W^ < a is satisfied, then the 
associated primary stationary solution of (2) is asymtotically stable. 

Proof. Let V® be the primary stationary solution of (2), as fiS' m Wlf < a is satisfied. Let also V M be the 
unique solution of the same equation with some initial condition V M (0) = 4> G J 7 , see theorem 3.2.1. We 
introduce a new function X = V^ — V® which satisfies: 

f d t X(z,A,t) = -aX(z,A,t) + f nxR+ W m (d 2 (z,z'),\\og(£)\)<d(X(z',A',t))f^dm(z') 
\ X(z,A,0) = cj ) (z,A)-V°(z,A) 
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where W m (d 2 (z,z'),\\og(^ 7 )\) = S' m W (d 2 (z,z'), |log(^)|) and the vector @(X(z,A,t)) is given by 
Q(X(z,A,t)) = S{^V^z,A,t)) - S(tiV°(z,A)) with S = (S' m )- l S. We note that, because of the 
definition of 6 and the mean value theorem \<d(X(z, A, t))\ < /i\X(z, A, t)\. This implies that |0(r)| < \r\ 
for all rei 



d t X(z,A,t) = -aX(z,A,t)+ f W m (d 2 (z,z / ), 



io g (^: 



=>dt [e at X(z,A,t) 



W m [<k{z,z% 



log(^) 



)&(X(z',A',t)) C ^dm(z') 
")e(X(z',A',t))^dm(z') 



X(z, A, t) = er at X(z, A, 0)+ f e -<*('-«) f W m (d 2 (z, z'), 



log(^) 



dA' 



Q(X(z', A', u))-T7 dm(z')du 



\X(z,A,t)\ <e~ ttt |A(z,A,0)|+/i / e - a ('-«) 



W m d 2 (*,<), 



!og(^) 



dA' 
|A-(z',A',u)|— ■ dm(z')du 



^||^(i)||oo<e- Qt ||^(0)||oo+^o 9 ^ / e-^-^HJC^IUdo 

Jo 

If we set: G(t) = e Q *||X(t)|| 00 , then we have: 

G(t) < G(0) + A*W 9 ^ n / G(u)du 
Jo 

and G is continuous for all t > 0. The Gronwall inequality implies that: 

G{t) < G(0)e» w ° s ™' 
=>||A-(*)|| 00 <e^ w o^-a)t|| A - (0) || oo) 
and the conclusion follows. 

5 Spatially localised bumps in the high gain limit 

In many models of working memory, transient stimuli are encoded by feature-selective persistent neural 
activity. Such stimuli are imagined to induce the formation of a spatially localised bump of persistent 
activity which coexists with a stable uniform state. As an example, Camperi and Wang [ ' I] have proposed 
and studied a network model of visuo-spatial working memory in prefontal cortex adapted from the ring 
model of orientation of Ben-Yishai and colleagues [3] . It is therefore natural to study the emergence of 
spatially localised bumps for the structure tensor model in a hypercolumn of VI. We only deal with the 
reduced case of equation (7) and keep the general case for future work. 

In order to construct exact bump solutions, we consider the high gain limit fx — > oo of the sigmoid 
function. As above we denote by H the Heaviside function defined by H(x) = 1 for x > and H(x) = 
otherwise. Equation (7) is rewritten as: 



□ 



d t V(z,t)=-aV(z,t)+ / W{z,z')H(V(z',t) - K)dm(z') + I(z,t) 



(20) 



= -aV{z,t)+ W(z,z')dm{z') + I{z) 

J{z'&a\v(z',t)>K,} 

We make the assumption that the system is spatially homogeneous that is, the external input I does 
not depend upon the variables t and the connectivity function depends only on the hyperbolic distance 
between two points of D: W(z, z') — W(d 2 (z, z')). We also introduce a threshold k to shift the zero of 
the Heaviside function. 
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5.1 Stationary pulses 

Our aim is to construct a hyperbolic radially symmetric stationary pulse. Let us first consider a general 
stationary pulse: 



aV(z)= W(z,z')dm(z')+I(z) 

J{z'€B\V(z')>k} 

We assume that the set K = {z <E ID|V(z) > k} C D is compact. We note M{z,K) the integral 
J K W(z, z')dm(z'). The relation V(z) = K holds for all z € dK. 

In order to calculate M, we use the Fourier transform. First we rewrite M as a convolution product: 

M(z,K)= f W(z,z')dm(z') = f W{z, z')l K {z')dm(z') = W * l K (z) 

where W(z) d = f W(d 2 (z,0)). 

In [17], Helgason proves an inversion formula for the hyperbolic Fourier transform and we apply this 

result to W. 

W(z) = -i- / W(A)$ A (z)Atanh(^A)dA (21) 

4tt J r 2 



Then, 



M(z,K) =W*l K (z) = — / W(A)$ A *l A -(»Atanli(-A)dA 

4tt Jr 2 



It appears that the study of M(z, K) consists in calculating the convolution product $> * \k(z). 
Let z = h ■ O for h € G we have: 



$ A * l K (z) = / 1* ( ff • 0)$A(g _1 • 2)rf.g 

= f l K (g-0)^ x (g- 1 h-0)dg 

Jg 



<G 

for all g,heG, <5>\{g- l h ■ O) = $\{hT l g ■ O) so that: 

$ A * l«r («) = / ljf (s • 0)*a(/i _1 5 ' 0)dg - / 1 x (z')$a(^ 1 • z')<Mz') 
jg jd 

K 

5.1.1 Study of M(z,K) when X = B h (0,w) 

We now consider the special case where K is a hyperbolic disk centered at the origin of hyperbolic radius 
u, noted S/j(0,w). 

First step We start by giving an explicit formula for a useful integral. 
Lemma 5.1.1. For all uj > the following formula holds: 

$\(z)dm(z) = 7rsinh(cj) 2 cosh(a;) 2 $l (L 1 ' 1) (w) 

B h (0,w) 

where Bh(0,uj) is the hyperbolic ball and: 

$ A Q '%) = F(i(p + iA), i(p - iA); a + 1; - sinh(^) 2 ), 
where a + /3 + 1 = p and F is t/ie hyper geometric function of first kind. 
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Proof. We write z in hyperbolic polar coordinates, z — tanh(r)e (see appendix D). We have: 

pUJ p2lT 

2 



/ $x(z)dm(z) = - J j <f>x(ta,nh{r)e M )smh(2r)drd9 

JB h (0,u) 



10 Jo 
Because of the above definition of $a, this reduces to 

■k I $A(tanh(r))sinh(2r)dr 
Jo 

In [17] Helgason proved that: 

$ A (tanh(r)) = F(v, 1 - v; 1; - sinh(r) 2 ) 
with v = |(1 + iX). We then use the formula obtained by Erdelyi in [12]: 

F(v, l-v;l;z) = ~ (zF{v, 1 - v; 2; z) 

Using some simple hyperbolic trigonometry formulae we obtain: 

sinh(2r)F(^, 1 - v; 1; - sinh(r) 2 ) = — ( sinh(r) 2 F(^, 1 - v; 2; - sinh(r) 2 ) J , 

from which we deduce 

/ $ \{z)dm(z) = 7rsinh(w) 2 F(^, 1 - v\ 2; - sinh(w) 2 ) 

JB h (0,u) 

Finally we use the equality shown in [ ] : 

F(a, b; c; z) = (1 - z) c - a - b F{c -a,c-b;c; z) 

In our case we have: a = v, 6=1 — i/, c = 2 and z = — sinh(aj) 2 , so 2 — ^ = |(3 — zA), l + v= |(3 + iA). 
We obtain 

/ $ A (»dTO(z) = 7T sinh(w) 2 cosh(o;) 2 F( -(3 - JA), -(3 + *A); 2; - sinh(w) 2 ) . 

JB h {0,u) V 2 2 / 

Since Hypergeometric functions are symmetric with respect to the first two variables: 

F(a, b; c; z) = F(b, a; c; z), 
we write 

f(\$ ~ *A), ~(3 + tA); 2; - sinhH 2 ) = fQ(3 + *A), *(3 - *A); 2; - sinh(u;) 2 J = ^(u), 

which yields the announced formula 

<&\{z)dm(z) = 7rsinh(w) 2 cosh(u;) 2 $ A 1 ' 1 '(cL;) 



Bh(0,w) 



n 



Second step In this paragraph, we show that the mapping z = h ■ O = tanh(r)e* e — > f K <&\(h 1 • 
z')dm(z') is a radial function, i.e. it depends only upon the variable r. We then show that J„ , Q , <$>\(a,- r - 

z')dm(z') = $ A (a r • O) f BhM e^+V<z'^>dm(z>). 

Proposition 5.1.1. If z = h ■ O and z is written tanh(r)e* e with r = d2(z,0) in hyperbolic polar 
coordinates the integral J B , Q . <J>a(^ _1 ■ z')dm(z') depends only upon the variable r. 
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Proof. If z = tanh(r)e , then z — rotga r • O and h 1 = a_ r rot_g. Similarly z' — Tots>a r ' • O. We can 

write 

§x{h~ l ■ z')dm(z') = - / $ A (a- r rote'-ea r / • O) smh(2r' )dr'd9' 

h(0,w) 2 Jo ^0 

= - / / $A(a-r-rot^,a r / ■ 0)smh(2r')dr'dip = / $ A (a-r • z')dm(z'), 

2 Jo JO JB h (0,uj) 

which, as announced, is only a function of r. □ 

We now give an explicit formula for the integral J B ,„ , $ A (a_ r ' z')dm(z'). We first recall a formula 
from [17]. 

Lemma 5.1.2. For all g € G the following equation holds: 

* X (g- K z)= I e (-<A+D<»-0.6> e (*A+l)<*,6> d6 

J9D 

Proo/. See [17]. D 

It follows immediately that for all z e D and reRwe have: 

$ A (a_ r -z)= / e (-*A+l)<«r-0,6> e (»A+l)<*,6> d6 

J3B 

We integrate this formula over the hyperbolic ball Bh(0,ui) which gives: 

<S> x (a_ r -z')dm(z') = f I f e(- iX+ V <a r- < b >e( iX+ V <z '> b> db)dm(z'), 

B h (o,w) JB h {a,u>) \Jaa J 

and we exchange the order of integration: 

e (-iX+i)<a T -0,b>( f e {lX+ ^ <z '' b> dm(z')\db. 

S \JB h (0,u) J 

We note that the integral J B , Q * e A +!)< z ' b> dm(z') does not depend upon the variable b — e 1 ^ . Indeed: 



e 



(i\+l)<z',b> l 



^^^^rr L/wwy^U 2 S inh(2x)d^ 



B„(o,*) ' 2i j \\tmh(x)e i6 - e^\ 2 

1 [" (** ( 1 - tanh(a:) 2 \ ^ .,,„,, f/1 

ill < sinh(2x)dxdd 



2i i V|tanh(a0e«(9-*>-l| 

I r f 2 * ( 1 - tanh(2;) 2 



2 Jo Jo Vltan 11 ^) 610 ' -!| 2 
and indeed the integral does not depend upon the variable 6 



sm 



i.(2x)da:d0 / , 



e (*A+l)<z\ b>dm(/) = / e (' iA+1 ) <Z '' 1> dTO(z'). 

Bh(0,w) JB h (0,oj) 



Finally, we can write: 



$ A (a_ r • z')dm(z') = e (-iX+i)<a r -o,b> db e (' A+1 ) <2 ' 1> dm(z') 

B h (o,u) Job JB h (o,uj) 

= $_ A (a r -0)/ e ( - lX+ ^ <z '' 1> dm{z') = <5> x {a r -0) f e (,A+1)<zU> rfm(z'), 

JB h {0,u) JB h {0,u) 

because $ A = $- A (as solutions of the same equation). 
This completes the proof that: 

I $x(a-r-z')dm(z') = ®x(ar-0) f e {iX+x)<z ' ' x> dm{z') (22) 

JB h (0,ui) JB h (0,u) 
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The main result At this point we have proved the following proposition. 

Proposition 5.1.2. If K = Bh(0,oj) and z = rotea r -0 £ K, M(z,K) is given by the following formula: 

M(z,B h (0,uj)) d = f M(r,iu) = ±- f W(A)$ A (a, • 0)* A (w)Atanh( ^-\)d\ (23) 

where 

* A («) = f ( e ( lA+1 ) <z '' 1> d™(z') (24) 



We are now in a position to obtain an analytic form for A4(r,uj). It is given in the following theorem. 
Theorem 5.1.1. For all (r,w) e R + x R+ : 

M{r,Lo) = - sinh(w) 2 cosh(w) 2 f W(A)$ A °' 0) (r)$ A M) (w)Atanh(^A)dA (25) 

Proof. We prove 



*aH = / <f>x(z)dm(z). 

JB h (0,Lo) 

Indeed, in hyperbolic polar coordinates, we have: 

* A («) = 4/ u / 2 "e( lA + 1 )<™ t '» a '- O ' 1> sinh(2r)dr^ 

= I Jo" I? e( lA+1 ) <a - ' e ^ 8> sinh(2r) dr d9 
= * Jo Sao e {lX+1)<ar °' b> db sinh(2r) dr 
= 7r J " $ A (a r • O) sinh(2r) dr 



On the other hand: 



This yields 



= n Jo Jo* $ a (flr • O) sinh(2r) dr 



*a(w) = / $ A (z)dm(z), 

JB h (0,u) 



and we use lemma (5.1.1) to establish (25). □ 

Discussion Let us point out that our result can be linked to the work of Folias and Bressloff in [15]. 
They constructed a two-dimensional pulse for a general, radially symmetric synaptic weight function. 
They obtain a similar formal representation of the integral of the connectivity function w over the disk 
B(0,a) centered at the origin O and of radius a. Using their notations, 

/>27r />a />oo 

M(a,r)= / toflr - r'\)r' dr' dQ = 2vra / w(p)J {rp)J 1 {ap) dp 

Jo Jo Jo 



where J v {x) is the Bessel function of the first kind and w is the real Fourier transform of w. In our case, 



instead of the Bessel function, we find <&v'' v ' (r) which is linked to the hypergeometric function of the 



first kind, as explained in lemma 5.1.1. 

We next adapt the results proved by Folias and Bressloff in [ ] to the hyperbolic case. 



19 




0.1 0.15 

width omega 



Figure 1: Plot of N(tu) defined in (27) as a function of the pulse width u> for several values of the input 
amplitude X and for a fixed input width a = 0.05. The connectivity function has been set to w(r) = e~ i> , 
with b = 0.2. 



5.1.2 A hyperbolic radially symmetric stationary-pulse 

We note V(r) a hyperbolic radially symmetric stationary-pulse solution of (20) where V depends only 
upon the variable r and is such that: 

V{r) > k, re [0,w[, 

V(W) = K, 

V(r) < K, r e]w,oo[, 

V(oo) = 0. 
aV(r) =M(r,u) +I(r) 



and 

Substituting into (20) yields: 



where M(r,ui) is defined in equation (25) and I(r) = Xe 2^ is a Gaussian input. 
The condition for the existence of a stationary pulse is given by: 



<!<.} 



where 



an = M{uj) + I(u) = N(u) 



M(lu) d = f M{u,u) = - sinh(cj) 2 cosh(a;) 2 f W(A)$f 0) (w)^ 1,1) (w)Atanh(-A)dA 
4 ./» 2 



(26) 



(27) 



The function N(lu) is plotted in figure 1 for a range of the input amplitude X. The horizontal dashed 
lines indicate different values of an, the points of intersection determine the existence of stationary pulse 
solutions. 

We now show that for a general monotonically decreasing weight function w, the function M(r,cu) 
is necessarily a monotonically decreasing function of r. This will ensure that the hyperbolic radially 
symmetric stationary-pulse solution (26) is also a monotonically decreasing function of r in the case of 
a Gaussian input. Differentiating M. with respect to r yields: 



dM, , 1 

^T (r ' w) = 2 



d_ 

Or 



We have to compute 
d 



Or 



«)((i2(tanh(r),tanh(r')e ie )) sinh(2r' )dr'd6 







w(d 2 (tanh(r),tanh(r')0) ) = w'(d 2 (tanh(r),tanh(r')e^))— ( d 2 (tanh(r), tanh(/)e ltf ) 
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It is result of elementary hyperbolic trigonometry that 

/ tanh(r) 2 + tanh(r') 2 — 2tanh(r) tanh(r') cos(6) 



d 2 (tanh(r),tanh(r')e ltf ) = tanh" 



1 + tanh(r) 2 tanh(r') 2 — 2tanh(r) tanh(r') cos(9) 

we let p = tanh(r), p' — tanh(r') and define 

p 2 + p' 2 - 2pp' cos(9) 



f p -Ap) 



1 + p 2 p 12 - 2pp' cos(6) 



It follows that 



a toh- (Ji^ - ^ M 



and 



dP * V V " *') *{l-FprM)V*7M' 

8_ 2(p-fJco S (O)) + 2pp'(f/-pco8(0)) 

dp p '' e{P) {l + p 2 p ,2 -2pp'cos(6)) 2 

We conclude that if p > tanh(w) then for all < p' < tanh(w) and < 9 < 2ir 

2{p-p'cos{9)) + 2pp'(p' -pcos{9)) > 0, 

which implies A4(r, oj) < for r > w, since w' < 0. 

To see that it is also negative for r < u, we differentiate equation (25) with respect to r: 

^(r,w) = \ sinh(w) 2 coshH 2 / W(A)^-^°' 0) (r)$^ 1) (w)Atanh(jA)dA 
or 4 jr c*r 2 

The following formula holds for the hypergeometric function (see Erdelyi in [12]): 

— F(a, 6; c; z) = — F(a + l,b + l;c + l;z) 
dz c 

It implies 

^-$f 0) (r) = -isinh(r)cosh(r)(l + A 2 )$^ 1) (r). 

Substituting in the previous equation giving ^£f- we find: 

^(r, W ) = -^ sinh(2^) 2 sinh(2r) | W(A)(1 + A 2 )$^ 1 )(r)$^ 1) HAtanh(|A)rfA, 

implying that: 

dM dM 

sgn(— — (r,w)) =sgn(—— (w,r)). 
or ar 

Consequently, ^^(r, w) < for 7' < w. Hence 1/ is monotonically decreasing in r for any monotonically 
decreasing synaptic weight function w. 

5.1.3 Linear stability analysis 

We now analyse the evolution of small time-dependent perturbations of the hyperbolic stationary-pulse 
solution through linear stability analysis. 

Spectral analysis of the linearized operator Equation (20) is linearized about the stationary 
solution V(r) by introducing the time-depndent perturbation: 

v(z,t) = V(r) + <P{z,t) 

This leads to the linear equation: 

d t <f>(z, t) = -a<f>{z, t)+ J W{z, z')H'{V{r') - k)<j>{z' , t)dm(z'). 
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We separate variables by setting cf>(z,t) = (j)(z)e l3t to obtain the equation: 

(J3 + a)(j)(z) = f W(z, z')H'(V(r') - n)<j)(z')dm(z') 



Introducing the hyperbolic polar coordinates z — tanh(r)e and using the result: 

H'(V(r) - «) = 8{V(r) - k) = fc-^ 

we obtain: 

((3 + a)(/)(z) = I f f VK(tanh(r)e ^ ^tanh(r0e ^e ')^f^(/)(tanh(r')e ^e ')sinh(2/)dr'd6'' 
2 Jo Jo \v Ml 

smh{2 " :] f W r (tanh(r)e je ,tanh(w)e ie ')0(tanh(w)e l9 ')d6l' 



With a slight abuse of notation we are led to study the solutions of the integral equation: 

((3 + a)cf>(r, 9) = f^ft| J * W(r, w; 0' - fl)^(w, 0')^' (28) 

where: 

/ tanh(r) 2 + tanh(w) 2 — 2 tanh(r) tanh(w) cos((/?) \ 



W(r, uj; ip) = w o tanh 



1 + tanh(r) 2 tanh(o;) 2 — 2 tanh(r) tanh(w) cos(ip) 



Essential spectrum If the function <p satisfies the condition 

r 2ir 

W(r,uj;e')(f>(u,9-e')de' = Vr, 

o 

then equation (28) reduces to: 

f3 + a = 

yielding the eigenvalue: 

/3 = -a < 

This part of the essential spectrum is negative and does not cause instability. 

Discrete spectrum If we are not in the previous case we have to study the solutions of the integral 
equation (28). 

This equation shows that <f>(r, 9) is completely determined by its values </>(w, 9) on the circle of equation 
r = ui. Hence, we need only to consider r — u>, yielding the integral equation: 

[fi + a)4>{u,, 9) = ^j^j I * W(w, uj; 9')cj>{uj, 9 - 9')d9' 
z \ v Ml Jo 

The solutions of this equation are exponential functions e 7 , where 7 satisfies: 

By the requirement that (j> is 27r-periodic in 9, it follows that 7 = in, where n € Z. Thus the integral 
operator with kernel W has a discrete spectrum given by: 

2tt 



^^^ww\!J m ^ el) ^ ne ' de ' 



sMl ^ r ffi tanh -i A 2tanhM 2 (l-cos( ? ) ^^ 



2|V(w)| J I V l+tanh(cj) 4 -2tanh(w) 2 cos(6") 
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sinh(2w) r & tanh -! 2^(^(6') e _ i2nS , de , 



|V'(w)| Jo \ V(l - tanh(w) 2 ) 2 + 4tanh(w) 2 sin(0') 

/?„ is real since: 

^,„ s sinh(2u;) f n , , / 2 tanh(w) sin(6»') \ , 

Qf(S„) = — , — ^-f / tD o tanh -1 i ; i ; — sin(2n(9' )dd' = 0. 

|F'(w)| i \ y/(l - tanh(w) 2 ) 2 + 4 tanh(w) 2 sin(0') 2 / 

Hence, 

_ /0 , sinh(2w) r _ ,_ x ( 2tanh(u;)sin(0 / ) \ ,„„,,„, 

P n = U(B n ) = -a + , T „; ,/ / w o tanh 1 \ > \ J = co8 (2n6)d6' 

\V'(u)\ J \V( l - tanh(u;) 2 ) 2 + 4 tanh(w) 2 sin(6") 2 / 

Since w o tanh - (r) is a positive function of r, it follows that: 

/3„</3 
Stability of the hyperbolic stationary pulse requires that for all n > 0, /3 n < 0. This can be rewritten as: 

sinh(2w) r _ , i / 2tanh(w)sin(0') \ . „., 7 „, 

-s — hvr / * ° tanh -1 v ' I ' cos(2n0')d6' <a n>0 

|V'(w)| 7 V V(l - tanh(w) 2 ) 2 + 4 tanh(w) 2 sin(fl') 2 / 

Using the fact that j3 n < fio for all n > 1, we obtain the reduced stability condition: 

WoM 
|V'( W )| a 

where 

,„, ,,*/,,,„>/". , _i / 2tanh(w)sin(<9') \ ,„. 

VVUw) = sinh(2w) / to o tanh x - y v ; \d6' 

Jo \ y/{\ - tanhH 2 ) 2 + 4 tanh(o;) 2 sin(6>') 2 j 

From (26) we have: 

V*(w) = ~(-M r (w) + J'(w)) 
a 

where 

Mr(u)^ -~(w,w) = ~sinh(2w) 3 /w(A)(l + A 2 )$^ 1 ' 1) (w)$^ 1) (w)Atanh(^A)dA 
or 64 J K ' 2 

We have previously established that A4 r (w) > and /'(w) is negative by definition. Hence, letting 
T>(u) = \I'(u)\, we have 

\V (oj)\ = -(M r {ui) + V{uj)). 
a 

By substitution we obtain another form of the reduced stability condition: 

V(uj) > Wo(lo) - M r {uj) (29) 

We also have: 

X (w) = — A4(w,w) = — (w,w) + -g— (w,w) = W (w) - M-M, 

and 

N'{u) = .M'(w) + J'(w) = Wo(w) - .Mr(w) - V{lu), 

showing that the stability condition (29) is satisfied when N'(uj) < and is not satisfied when N'(ui) > 0. 
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6 Numerical results 

The aim of this section is to numerically solve (7) for different values of the parameters. This implies 
developing a numerical scheme that approaches the solution of our equation, and proving that this scheme 
effectively converges to the solution. 

Since equation (7) is defined on D, computing the solutions on the whole hyperbolic disk has same the 
complexity as computing the solutions of usual Euclidean neural field equations defined on R 2 . As most 
authors in the Euclidean case [15, 22, 23, 25], we reduce the domain of integration to a compact region 
of the hyperbolic disk. Practically, we work on the Euclidean ball of radius a — 0.5 and center 0. Note 
that a Euclidean ball centered at the origin is also a centered hyperbolic ball, their radii being different. 
We have divided this section in four parts. The first part is dedicated to the study of the discretization 
scheme of equation (7). In the following three parts, we study the solutions of different connectivity 
functions: exponential function 6.2, Gabor function 6.3 and a difference of Gaussians function ??. 

6.1 Numerical schemes 

Let us consider the modified equation of (7): 

d t V(z,t) = -aV(z,t) + J B{0a) W(z,z')S(V(z',t))dm(z') + I(z,t) t£j 

V(z,0) = V (z) ( ' 

We assume that the connectivity function satisfies the conditions (C1)-(C2). Moreover we express z in 
(Euclidean) polar coordinates such that z = re 10 , V(z, t) = V(r, 0,t) and W(z,z') = W(r, 9, r',6'). The 
integral in equation (30) is then: 

r'dr'de' 
W(z,z')S(V(z',t))dm(z')= I / W(r, 0, r' , 6')S(V{r', 0' , t))- / 

B(0,a 



W(z,z')S(V(z',t))dm(z')= / W(r,e,r',6')S(V(r',6',t))- n , , ., 

Jo Jo I 1 — r ) 



We define K to be the rectangle K d = [0, a] x [0, 2tt]. 



6.1.1 Discretization scheme 



We discretize 1Z in order to turn (30) into a finite number of equations. For this purpose we introduce 
hi = f , N £ N* = N\{0} and h 2 =|,Msff*, 

V»e[i,JV + i] ri = {i-i)hi, 

Vje[l,M + l] 9 3 = (j-l)h 2 , 

and obtain the (N + 1)(M + 1) equations: 

dV f r'dr'dO' 

^(r l ,e j ,t) = -aV(n,9 J ,t) + j^W(n,9 J ,r / ,e / )S(V(r / ,e\t))^—-^ +1(^9^1) 

which define the discretization of (30): 



I <^(t) = -aV(t)+W-S(V)(t)+I(t) teJ 
\ V(0) = V 
where V(t) E A / Jat+i.m+i(R) 1 , V(t)i,j — V(ri,9j,t). Similar definitions apply to I and Vq. Moreover: 



(31) 



w • s(v)(t) itj = J w( n , Ojj, e')s(v(r', 9', *)) (1 _ r , 2)2 

It remains to discretize the integral term. For this as in [14], we use the rectangular rule for the quadrature 
so that for all (r, 9) £ 1Z we have: 



r a r 2n r'dr'dO' N+1M+1 

/ / w(r,ey,d')s(v(r',e',t))^ — ^2=^2 V V w^e^eosiv^e^t))-^.,. 

Jo Jo U-r ; k=1 l=1 V-r k ) 



2\2 



Mn,p(S-) is the space of the matrices of size n X p with real coefficients. 
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We end up with the following numerical scheme, where Vij(t) (resp. 2kj(t)) is an approximation of 
Vij(t) (resp. Iij), V(t,j) € [l.JV + 1] x [1.M + 1]: 

, JV+1M+1 ~ 

fc=i j=i 
with WJg =/ T^(r„ ^,r fe , 0j)pz%5*- 

6.1.2 Discussion 

We discuss the error induced by the rectangular rule for the quadrature. Let / be a function which is 

C 2 on a rectangular domain [a, b] x [c,d]. If we denote by Ef this error, then \Ef\ < - — 4 „ ln ||/||c 2 

where m and n are the number of subintervals used and ||/||c 2 = E| a |<2 su Pk6]x[c,(J \d a f\ where, as 

usual, a is a multi-index. As a consequence, if we want to control the error, we have to impose that the 

solution is, at least, C 2 in space. 

Four our numerical experiments we use the specific function ode45 of Matlab which is based on an explicit 

Runge-Kutta (4,5) formula (see [ '] for more details on Rung-Kutta methods). 

We can also establish a proof of the convergence of the numerical scheme which is exactly the same as in 

[14] excepted that we use the theorem of continuous dependence of the solution for ordinary differential 

equations. 

6.2 Purely excitatory exponential connectivity function 

In this subsection, we give some numerical solutions of (7) in the case where the connectivity function 
is an exponential function, w(x) — e~T~, with b a positive parameter. Only excitation is present in this 
case. In all the experiments we set a = 0.1 and S(x) = 1+ l-^ with /i = 10. 

Constant input We fix the external input I(z) to be of the form: 

d 2 (z,0) 2 
I{z) = le v 2 

In all experiments we set I = 0.1 and a — 0.05, this means that the input has a sharp profile centered 
at 0. We show in figure 2 plots of the solution at time T — 2500 for three different values of the width b 
of the exponential function. When 6=1, the whole network is hightly excited, see figure 2(a). When b 
changes from 1 to 0.1 the amplitude of the solution decreases, and the area of high excitation becomes 
concentrated around the external input. 

Variable input In this paragraph, we allow the external current to depend upon the time variable. 

We have: 

_ d 2 (z,Zo(t)) 2 

I(z, t) = le~ a 2 

where Zo(^) = roe tnot . This is a bump rotating with angular velocity Qq around the circle of radius Tq 
centered at the origin. In our numerical experiments we set r = 0.4, J7 = 0.01, X = 0.1 and a = 0.05. 
We plot in figure 3 the solution at different times T = 100, 150, 200, 250. 

High gain limit We consider the high gain limit [i — > oo of the sigmoid function and we propose to 
illustrate section 5 with a numerical simulation. We set a = 1, K = 0.04, u — 0.18. We fix the input to 
be of the form: 

d 2 (z,0) 2 

I(z) = le <J 2 

with X — 0.04 and a — 0.05. Then the condition of existence of a stationary pulse (27) is satisfied, see 
figure 1. We plot a bump solution according to (27) in figure 4. 
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(c) 6 = 0.1 



Figure 2: Plots of the solution of equation (7) at T = 2500 for the values /i = 10, a = 0.1 and for 
decreasing values of the width b of the connectivity, see text. 
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Figure 3: Plots of the solution of equation (7) in the case of an exponential connectivity function with 
b = 0.1 at different times with a time-dependent input, see text. 
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(a) b = 0.2 



(b) b = 0.2 



Figure 4: Plots of a bump solution of equation (26) for the values a = 1, K = 0.04, u = 0.18 and for 
b = 0.2 for the width of the connectivity, see text. 



6.3 Excitatory and inhibitory connectivity function 

We give some numerical solutions of (7) in the case where the connectivity function is a Gabor function. 
In all the experiments we set a = 0.1 and S(x) = l+ l~^ with /i = 10. The Gabor function is given by 

G(x) = -4= 1 1 — 2w ) e~~s~, where b is a positive parameter, takes positive values near the origin and 
negative values further away. 

The external input I(z,t) is equal to zero. In figure 5, we show several solutions of equation (7) 
at T = 2500 for decreasing values of the width b of the connectivity. We see the emergence of multi- 
bump solutions. In figure 5(a), there is one bump centered at O. When decreasing the width b of the 
connectivity, for the range of values b £ [0.2; 0.35], the solutions present an increasing number of areas 
of high level of activity. For example, in figure 5(d), there exist five such areas. Figure 5(c) shows the 
appearance of a second crown of localized patterns centered at 0, and figure 5(f) shows three such crowns. 
Note that all the solutions display interesting symmetries. 

We slightly change the shape of the connectivity function by choosing the difference of two Gaussians 

f e "a . In all the experiments, we take U\ = 0.1, <r 2 = 0.2 and A = 1. 



As shown in figure 6 the two connectivity functions, Gabor and difference of Gaussians, feature an 
excitatory center and an inibitory surround. We illustrate the behaviour of the solutions when increasing 
the slope \i of the sigmoid. We set the sigmoid S(x) = 1+ "L^ — \ so that it is equal to at the origin 
and we choose the external input equal to zero, I(z, t) = 0. In this case the constant function equal to 
is a solution of (7). 

For small values of the slope fi, the dynamics of the solution is trivial: every solution asymptotically 
converges to the null solution, as shown in figure 7(a). When increasing /i, the stability bound, found in 
subsection 4.5 is no longer satisfied and the null solution may no longer be stable. In effect this solution 
may bifurcate to other, more interesting solutions. We plot in figures 7(b), 7(c), 7(d), 7(e) and 7(f), some 
solutions at T = 2500 for different values of /i. We can see exotic patterns which feature some interesting 
symmetries. The formal study of these bifurcated solutions is left for future work. 



7 Conclusion 

We studied the existence and uniqueness of a solution of the equation of a smooth neural field model of 
structure tensors to model the representation and processing of texture and edges in the visual area VI. 
We also detailed the analysis of the stationary solutions of this nonlinear integro-differential equation. 
In both cases functional analysis and the theory of ordinary differential equations have allowed us to 
introduce the framework in which our equations are well-posed and to begin to characterize their solu- 
tions. This is interesting and important in itself but should also be useful for future investigations such 
as bifurcation analysis. 
We have completed our study by constructing and analysing spatially localised bumps in the high-gain 
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Figure 5: Plots of the solutions of equation (7) in the case of a Gabor connectivity function at time 
T = 2500 for the values fj, = 10, a = 0.1 and for decreasing values of the width b of the connectivity, see 
text. 
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Figure 6: Plot of the Gabor function G with b = 0.2 and in dashed lines the difference of Gaussian 
functions w with o\ = 0.1, oi = 0.2 and A = 1, see text. They both feature an excitatory center and an 
inhibitory surround. 

limit of the sigmoid function. It is true that networks with Hcavisidc nonlincarities are not very realistic 
from the neurobiological perspective and lead to difficult mathematical considerations. However, taking 
the high-gain limit is instructive since it allows the explicit construction of stationary solutions which is 
impossible with sigmoidal nonlinearities. We constructed what we called a hyperbolic radially symmetric 
stationary-pulse and presented a linear stability analysis adapted from [15]. 

Finally, we illustrated our theoretical results with numerical simulations based on rigorously defined 
numerical schemes. We hope that our numerical experiments will lead to new and exciting investiga- 
tions such as a thorough study of the bifurcations of the solutions of our equations with respect to such 
parameters as the slope of the sigmoid and the width of the connectivity function. 
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A Volume element in structure tensor space 

Let T be a structure tensor 



T 



A 2 its determinant, A > 0. T can be written 



Xi x 3 
x 3 x 2 



T=AT, 



where T has determinant 1. Let z = Z\ + iz 2 be the complex number representation of T in the Poincare 
disk D. In this part of the appendix, we present a simple form for the volume element in full structure 
tensor space, when parametrized as (A,z). 



Proposition A. 0.1. The volume element in (A, z\ 1 Z2) coordinates is 



dV = 



dA dz\ dzi 



(32) 



Proof. In order to compute the volume element in (A,Z!,z 2 ) space, we need to express the metric gj- 
in these coordinates. This is obtained from the inner product in the tangent space Tj- at point T of 
SDP(2). The tangent space is the set S(2) of symmetric matrices and the inner product is defined by: 

g T (A,B)=tr(T- 1 AT- 1 B), A,BeS(2), 
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(a) fj, = 1 



(b) M : 





(C) (l: 



(d) fj. = 10 





(f) /j = 30 



Figure 7: Plots of the solutions of equation (7) in the case where the connectivity function is the difference 
of two Gaussians at time T — 2500 for a — 0.1 and for increasing values of the slope /x of the sigmoid, 
see text. 
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We note that gr(A, B) = gf(A, B)/A 2 . We note g instead of gf. A basis of Tj- (or Tf for that matter) 
is given by: 

d __ r i o 1 _9__rooi o _r o i ~ 

and the metric is given by: 

5 -" = ^ ( ^'^- ) ' M ' = 1 ' 2 ' 3 

The determinant G-y of 517- is equal to G/A 6 , where G is the determinant of g = gf. G is found to be 
equal to 2. The volume element is thus: 

dV = -r^r dx\ dx2 dx3 
A 6 

We then use the relations: 

x\ = Aii X2 — Ax2 X3 = AX3, 

where ii, i = 1,2,3 is given by: 

x ' 2 = X-zt-zi 

ggg 
1-^1-^2 

The determinant of the Jacobian of the transformation (x±, X2, X3) — > (A, zi, 22) is found to be equal to: 

8A 2 
~(1-M 2 ) 2 
Hence, the volume element in (A, 21,22) coordinates is 



~~ \-z\-z\ 
X 3 - T^i^iJ 



n 



B Global existence of solutions 

Theorem B.0.1. Let O be an open connected set of a real Banach space JF and J be an open interval 
o/R. We consider the initial value problem: 



V'(t) = f(t,V(t)), 
V(t ) = V 



(33) 



FFe suppose that f G C (J x 0, J 7 ) and is locally Lipschitz with respect to its second argument. Then for 
all (to, Vo) E J x O, there exists r > and V G C 1 (]t r , to + t[, 0) unique solution of (33). 

Lemma B.0.1. Under hypotheses of theorem B.0.1, if Vi G C 1 (Ji,C) and V 2 G C 1 (J2,C) are two 
solutions and if there exists t G J\ PI J2 suc/i i/iat Vi(<o) = V 2 (£o) i/ien: 

Vi(t) = V 2 (t) /or a// 1 G Ji n J 2 

This lemma shows the existence of a larger interval Jo on which the initial value problem (33) has a 
unique solution. This solution is called the maximal solution. 

Theorem B.0.2. Under hypotheses of theorem B.0.1, let V G C (Jq,0) be a maximal solution. We 
note by b the upper bound of J and j3 the upper bound of Jq . Then either )3 — b or for all compact set 
K C O, there exists n < j3 such that: 

V(t) G 0\JC, for allt>r/ with t G J a - 

We have the same result with the lower bounds. 

Theorem B.0.3. We suppose f G C (J x J 7 , J 7 ) and is globally Lipschitz with respect to its second 
argument. Then for all (to, Vo) G J x T , there exists a unique V G C l (J,F) solution of (33). 
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C Proof of lemma 3.1.2 

In this section we prove the following lemma. 

Lemma CO. 2. When W is a mexican hat function of d 2 (T ,T') and independent oft, then: 

3 

W=7 Y ( CT i e2ff W(\/2^i) -Ao- 2 e 2 °hrf(y2a 2 ^ 
where erf is the error function defined as: 

2 f x _ 2 
erftx) = — = / e u du 

V 71 " Jo 

Proof. We consider the following double integrals: 

(log A -log A') 2 d 2 {z,z') 
+ °° r 1 " ^ - 9rr 2 dA' dz[dz'o , , 



V^^ 2 A'(l-|*'l 2 ) 2 

so that: 

TF = Hi - AE 2 

Since the variables are separable, we have: 

v'\2 All-, J\ 



:f 



(logA-logA') 2 , , di( 2 ,z') 



+ °° 1 "" o 72 dA '\( f 2cr 2 dz i dz 2 



/o v^of A '/V' /d {l-\A 2 Y 

One can easily see that: 



(log A -log A') 2 
+°° 1 ~2 dA' _ 1 

r-i ~~ 2rr 2 

We now give a simplified expression for Sj. We set /i(a;) = e t and then we have, because of 

lemma 3.1.1: 

S, - 4h / /,(d 2 (0,z'))dm(z') = i / /,(arctanh(|z'|))dm(z') 

V 2 JD V 2 Jd 
1 f 1 f 2lT rdrdO r- f 1 rdr 



/i(arctanh(r))— — - = VZtt / /i(arctanh(r)) 



VlJo Jo ' ' " "{l-r*Y J J ' v v "{1-r^Y 

arctanh (r) 
= V2n I e 



2cr 2 rdr 







(l_ r 2)S 



dr 

The change of variable x — arctanh(r) implies dx = and yields: 

1 — r 1 



Sj = v^tt / e 2cr ? 2 — dx = V2tt / e 2<j2 sinh(a;) cosh(a;)dj: 

Jo 1 — tanh (x) Jo 



2 x 2 x 2 

+oo — — — =■ / /-+oo — — — ^-+2k r+oo — — — k — 2x 



X 



-^ / e 2fJ * sinh(2a:)da: = — ^ / e 2ct * efa - / e 2ct * dx 



2^2 e '^io 



(*-2a 2 ) 2 (* + 2a 2 ) 2 



+oo — — = /■+oo 



e 2cr2 dx- / e 2ff2 da; 
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□ 



then we have a simplified expression for H$: 

3 

H« = ^-(Tie 2 ^ erf (V2ai) 

D Isometries of D 

We briefly descrbies the isometries of D, i.e the transformations that preserve the distance di- We refer 
to the classical textbooks in hyperbolic goemetry for details, e.g, [20]. The direct isometries (preserving 
the orientation) in D are the elements of the special unitary group, noted SU(1, 1), of 2 x 2 Hermitian 
matrices with determinant equal to 1. Given: 

7 = ( | I ) such that H 2 - |/3| 2 = 1, 

an element of SU(1, 1), the corresponding isometry 7 in D is defined by: 

1 . z = 4, zeD (35) 

pz + a 

Orientation reversing isometries of D are obtained by composing any transformation (35) with the re- 
flexion k : z — > z. The full symmetry group of the Poincare disc is therefore: 

U(l,l) = SU(1,1)Uk-SU(1,1) 

Let us now describe the different kinds of direct isometries acting in D. We first define the following one 
parameter subgroups of SU(1, 1): 

K *{»**=( f ° e _ 4 V^eS 1 } 

.defr ( coshr sinhr \ , 

^ sinhr coshr J ' •* 

,Je/ r 1 1 + is —is \ _.-, 

N = ' {n s = . . ,seR 

V is I — is J ' 

Note that rot^ ■ z — e l ^z and also a r ■ O — tanhr. 

The group K is the orthogonal group 0(2). Its orbits are concentric circles. It is possible to express 

each point z £ D in hyperbolic polar coordinates: z = rot^a,. • O = t&nhre 1 ^ and r = d2(z,0). 

The orbits of A converge to the same limit points of the unit circle <9B, b± = ± when r — ¥ ±00. They 

are circular arcs in D going through the points 61 and b—\. 

The orbits of N are the circles inside D and tangent to the unit circle at b\. These circles are called 

horocycles with base point b±. N is called the horocyclic group. It is also possible to express each point 

z € ID) in horocyclic coordinates: z — n s a r ■ O, where n s are the transformations associated with the 

group N (s £ K) and a r the transformations associated with the subroup A (r £ K). 

Iwasawa decomposition The following decomposition holds, see [19]: 

SU(1,1) = KAN 

This theorem allows us to decompose any isometry of D as the product of at most thee elements in the 
groups, K, A and N. 
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